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ABSTRACT 


The  stochastic  nature  of  the  target  return  depends  on  numerous  factors: 
the  spatial  randomness  of  the  forest?s  EM  parameters,  their  temporal 
fluctuations ,  the  target  size  and  its  state  of  motion.  Our  attention  is  focused 
upon  the  calculation  of  the  mean  power  and  the  temporal  spectrum  of  the 
target  return,  quantities  of  direct  and  obvious  relevance.  The  spectral 
modifications  are  traceable  to  two  first  order  factors:  (  1  )  the  motion  of  the 
scattering  centers  within  the  random  vegetation  slab  and  (2)  the  target  motion 
through  a  randomly  inhomogeneous  field,  with  the  latter  likely  to  dominate. 

The  analysis  is  carried  out  within  the  framework  of  the  so-called 
’’distorted  wave  Born  approximation.  n  The  incoherent  (random)  scattering 
is  accounted  for  to  the  order  of  single  scatter,  while  coherent  (background) 
effects  such  as  refraction  at  the  air -vegetation  interface  and  ground  reflec¬ 
tions  are  properly  considered.  The  rather  laborious  analytical  sequence 
leads  (via  a  set  of  carefully  listed  simplifying  assumptions)  to  final  results 
which  are  analytically  simple,  numerically  tractable  and  their  physical  con¬ 
tent  readily  interpretable.  Two  major  assumptions  concerning  target  char¬ 
acteristics  should  be  mentioned:  (1)  the  study  dealt  exclusively  with  point 
targets  and  (2)  the  target  is  presumed  to  move  toward  or  away  from  the 
radar  with  constant  velocity. 
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A.  INTRODUCTION 


The  stochastic  nature  of  the  target  return  depends  on  numerous  factors, 
the  spatial  randomness  of  the  EM  parameters  and  their  temporal  fluctuations, 
the  properties  of  existing  interface  boundaries  and  those  of  the  target  and  its 
state  of  motion.  While  virtually  all  wave  statistics  are  relevant  under  appro¬ 
priate  circumstances,  our  attention  here  is  focused  upon  the  temporal  spec¬ 
trum  of  the  target  return  ,  a  quantity  of  direct  interest  for  design  and  opera¬ 
tional  considerations  of  radar  systems  with  doppler  capabilities.  The  spec¬ 
tral  modifications  of  a  narrow  band  signal  are  traceable  to  two  first  order 
factors:  (a)  the  time  dependence  characterizing  the  fluctuating  medium  (the 
vegetation  slab  in  our  application)  and  (b)  the  motion  of  the  target  through  a 
randomly  inhomogeneous  field,  with  the  latter  likely  to  dominate. 

The  objective  of  this  analytical  work  is  to  relate  the  desired  spectral 
information  to  the  presumably  known  statistical  measures  which  characterize 
the  random  slab  (specifically  its  space-time  correlation  function  or  the  cor¬ 
responding  Fourier  transform).  Several  avenues  are  open  to  us  within  the 
framework  of  the  previously  introduced  continuum  model'  ^  Of  these,  the 
simplest  alternative,  that  involving  a  single  random  scatter  event  of  either 
the  incident  or  the  return  wave,  is  to  be  pursued.  The  included  effects  are 
depicted  symbolically  in  Fig.  1,  in  which  refractive  and  ground  effects  are 
shown  and  are  indeed  accounted  for  in  the  analytical  sequence  that  follows. 

The  methodology  resembles  that  presented  in  Ref.  (1).  Conclusions  and 

The  analogous  problem  concerning  the  spectrum  and  other  statistical 
properties  associated  with  the  clutter  return  have  been  investigated  in  some 
detail.  The  results  and  conclusions  have  been  reported  in  Ref.  1. 
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Fig.  1.  Ray  configurations. 
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results  presented  in  it  are  not  repeated'  . 

The  conceptual  simplicity  inherent  in  the  (distorted  wave)  single  scatter 
theory  to  be  pur  sued  does  not  fully  propagate  through  the  analytical  sequence  which 
is  rather  laborious.  However,  through  the  imposition  of  various  simplfying 
assumptions,  most  of  which  not  too  objectionable,  the  final  result  (Eqs.  (68 
and  76))  is  analytically  simple,  numerically  useful  and  its  physical  content 
readily  inte rpr etable .  In  fact,  a  crude  but  direct  and  simple  heuristic  con¬ 
struction  of  the  final  result  is  presented  in  Sec.  F.  It  should  be  mentioned 
that  the  presented  analysis  retains  its  validity  and  applicability  beyond  the 
physical  context  of  immediate  interest. 

The  problem  at  hand  is  described  and  rigorously  formulated  in  Sec.  B 
where  a  formal  expression  (14)  for  the  target  generated  electric  field  (E^) 
at  the  receiving  antenna  is  derived.  An  explicit  expression  for  the  stochastic 
target  return  is  obtained  in  Sec.  C  consistently  with  thendis  tor  ted  wave11  Born 
approximation  in  which  background  refraction,  ground  effects  and  single  in¬ 
coherent  scattering  are  properly  accounted  for.  The  temporal  correlation 
coefficient  (C^fT))  and  the  associated  spectral  density  are  calculated  in  Sec. 

D.  The  main  result  (68)  is  obtained  subject  to  the  omission  of  the  motion 
associated  with  the  randomly  distributed  scattering  centers.  Sec.  D  termin¬ 
ates  with  several  interpreting  notes.  The  spectral  perturbation  associated 
with  the  random  motion  of  the  scattering  centers  is  presented  in  Sec.  E.  The 
main  result  is  contained  in  Eq.  (76).  The  content  of  Sec.  F  has  been  mentioned 

Reference  to  equations  in  Ref.  (1)  are  made  by  adding  I  (  )  in  front  of  the 

equation's  number. 
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above.  The  mean  power  return  is  calculated  in  Sec.  G  and  the  non-trivial 
question  regarding  the  definition  of  the  effective  scattering  volume  is  the 
subject  matter  of  Appendix  A. 

B.  FORMULATION  AND  DESCRIPTION  OF  THE  PROBLEM 

The  physical  configuration  comprises  a  "point11  target  moving  at  a 
constant  velocity  (v)  through  a  randomly  fluctuating  slab  region  depicted  in 
Fig.  1.  The  slab,  modeling  the  vegetation,  is  parameterized  electromagnet- 
ically  by  random  permittivity  fluctuations  (G  (r,  t))  superimposed  on  a  lossy 
background  represented  by  the  complex  permittivity  e7,  consistently  with 
the  "continuum"  model  presented  in  Ref.  (1).  A  heuristic  sequence  arguing 
the  validity  of  the  continuum  model  is  also  given  in  Ref.  (1)  and  need  not  be 
repeated.  The  target  is  characterized  by  an  effective  point  current  (dipole) 
element 


JT(r)  =  jo)  P^  -  0  E,  (r)  6  (r  -  r_(t)  ) 


(1) 


where  rr^(t)  denotes  the  target  location,  and  (3/jo)^  is  the  target  polarizability 


4tt 


0  e  -  1 
i  r 


(e.  g.  (3  =  -  — — r —  (k~  a)  - ^  for  a  small  dielectric  sphere  of  radius 

Li  £  £42 

o  2  r 

a  <<\  and  relative  permittivity  e  k^  the  wave  number  in  region  2).  The 
electric  fields  E^  (r,  t)  and  E^  (r,  t)  in  regions  1  and  2  respectively  are 
solutions  of  the  vector  wave  equations. 


2 

V  x  V  x  JS  -  kj  Jq(£)>  z  >  0  (2) 

2  2 

V  xVx  E~  -  k0  E~  -  uu  a  e(r,  t)  E0  =-ja>  u  J-,  z  <  0  (3) 

— z  z  — z  o  o  —  — z  o  o  —  i. 

together  with  the  continuity  of  the  transverse  electric  and  magnetic  fields  at 
z  =  0,  the  vanishing  of  the  transverse  electric  field  at  z  =  -h  (h  denoting  the 
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effective  slab  thickness)  and  the  radiation  condition  at  z  -»  ®  in  (2) 

denotes  the  exciting  current  density  situated  in  region  1.  With  the  Green's 
functions  (G^  r,  _r ')  )  defined  by 


Tx  -  k^G,  =  0.  z  >  0  (4) 

V  X  V  X  G-  -  k,2  cL  -  0)  2  n  0  (r,  t)  G  ,  =  I  6  (r  -  r  '  ),  0  <  z  <  h  (5) 

Eq s .  (2,  3)  may  be  formally  inverted.  One  obtains, 

Ej  T  (£  |  t)  =-ju)Q  Uq  J'd 3  rl  Gj  (r,  •  JT  ( r  j ) 

2  2 

=  -j«"o  Uo  5  Si  (1-  -T  (t)  I  ^  '  — 2^-T^'  I 


where  E,^  is  the  target  generated  field  defined  by 


— 1,  2T  —  1,  2  "  —  1,  2  {-T  0) 


(7) 


A 


The  difficulty  in  (6)  stems  from  the  fact  that  neither  G  ^  ?  nor  E?  are  known. 

The  ultimate  validity  of  the  theory  depends  on  the  accuracy  with  which  these 
quantities  can  be  estimated. 

A 

We  now  focus  attention  on  G]  Let, 

^  1  >  ^ 

G  L  (r,  r_j  |  t)  =  G  x  (_r,  )  +  g  j  (r_,  |  t) 

2  2  2 
A  (i  .7 

where  G  =  G  (G=  0)  and  g  is  calculated  via  the  distorted  wave  Born 
~  i  ^  ^  i  / 

2  2  2 

approximation.  We  have  (exactly) 


(8) 


Vx  Vx  £i  -  ki  £i  =  0 

y  x  v  x  g2  -  k  2  g2  =  o)o2  po  c  g2 


(9) 

(10) 
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(11) 


which,  upon  inversion,  results  in  the  integral  equation 

gj  (£.’  £'11)  =  %  J  Si  (L>  Iq)  *  S2  (Il’l'  I  fc) 

2  V  2 

where  a  detailed  definition  of  the  effective  scattering  volume  (V)  is  presented 

in  Appendix  A.  Consistently,  with  the  distorted  wave  Born  approximation,  the 
A 

unknown  in  (11)  is  replaced  by  G  which  is,  in  turn,  determined  from  the 
deterministic  background  problem^.  The  field  appearing  in  (6)  may  be 
written  as 


(o) 


E2  (r_  I  t)  =  Ez  (£)  +  E2s  (r_  I  t) 


(12) 


(o) 


where  E2  =  E2  (e=  0)  and 

—2s  ^lfc)  =  \  ^o  S  rl  C(-l’t}  S  <1-  l!>  •  S2  (rj  |t) 


(13) 


Once  again,  consistently  with  the  distorted  wave  Born  approximation,  E, 
in  (13)  is  to  be  replaced  by  E^°^.  The  substitution  of  Eqs.  (8,  12)  into  (6) 
results  in 

E1T  (1  |t)  =  -j»0  3  [Gj  (£,  £t)  -  E^0)  (£t)  +  g1  (£,  £t  |  t)  •  (£t) 

2  2  2 


+  Si  (£•  Lt)  •  52s  (rT  |  t)  ] 
2 


(14) 


where  each  of  the  three  terms  is  readily  identified  with  the  distinct  contri¬ 
butions  depicted  in  Figs,  la  and  lb,  and  the  double  scatter  term  has  been 
omitted. 

C.  THE  TARGET  RETURN 


Eq.  (14)  may  be  cast  into  the  form 

(o) 


E  ^  p  (£_  I  t)  “  E^  p  (£  I  t)  +  A  E.  p(_r  1 1) 

2  2  i 


(15) 
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where 


Ej°' (r_  |  t)  =  -juuo  0  Gj  (r_,  £t)  •  E^0’  (r_ 


(o) 


T) 


(16) 


and 


4E1T(r|t)=  -j%3  uo2  6  jd’rjClr,.,)  [G^r.r,).  O,  E<°’(rT 


(o). 


+  Gj  (r_,  £t)  -  G2  (£t>  £j)  -  E2°'  (r  j  )  j 

2 

,(o) 


(17) 


It  is  observed  that  the  term  corresponds  to  a  return  from  a  point 

2  ...  * 
target  moving  at  a  constant  velocity  through  a  uniformly  illuminated  region 

It,  therefore,  represents  a  spectrally  unperturbed  (with  the  exception  of  a 

simple  shift)  return.  E^0^  an<^  G  j  are  taken  to  be  these  derived  in  Ref.  (1) 

(Eqs.  1(98)  and  1(101),  respectively) 


z  [l  -  p  p  ]  2 

Sl(-’  -1}  ~  T~~  L—  2~  exP  {  -jkl  tl£  ■£!  I  +  7 

2rr  vlp  -  1  Lj  1 


|2."£i  I  ^Ti2.  i 


■  J  ]  +  cc- 


LLf2 

(18) 


—2°  *  (— 1 )  ~  ® 


2z  F(cp) 

°  '/t7Tli2 


exp  {-jkJU-  I  +  I  ^  ^  .,h 


(19) 

The  dyadic  Green's  function  contains  direct  ray  contributions  as  well  as 
contributions  associated  with  interface  (ground  vegetation  and  air  vegetation) 
reflections.  For  CL  air  vegetation  interface  reflections  may  be  safely  ignored 


The  transverse  field  inhomogeneity  associated  with  the  antenna  beam  may 
be  disregarded. 
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for  two  reasons.  Firstly,  the  interface  discontinuity  is  small  and  with  the 
exception  of  critical  incidence  so  is  the  associated  reflection  coefficient. 
Secondly,  since  the  medium  is  lossy  and  since  the  anticipated  targets  are 
situated  near  ground,  the  reflected  contribution  is  further  attenuated  due 
to  the  absorption  and  scattering  losses  and  become  insignificant  even  for 
larger  interface  discontinuity.  With  ground  effects  included,  one  has 

exp  [  -jk  |r  -£x|  -  a  |£  -  r '|  ]  / 

G?(r,  r)~ - ^ - (1  -  r  r  ),  r  =f'4r.t-  (20) 

~2 -  4  IT  |r  -  r'|  |  r-r’l 


and 


1/2 


exp{(-jk?  -  Ct)  [  |_p -_p '|  +  (z  +  z'+2h)  ]  } 


4Tr[|£-£/|  +  (z+z'+2h)  ] 


2,1/2 


^ol 


(21) 


where  rn  denotes  the  image-point  location  (Fig.  2).  Ground  effects  are 

accounted  for  by  replacing  G  and  in  Eq.  (17)  by  G,  +  G1t)and 

2  2 

E^°^  +  respectively.  The  subscript  R  denotes  ground  reflections 

which  for  the  horizontally  polarized  constituents  are  obtained  by  a  phase 
reversal  and  the  replacement  of  -z^  and  -z^  by  2h  +  z^  and  2h  +  z^, 
respectively.  It  has  been  shown  in  Ref.  (1)  that  single  ground  reflections 
are  properly  accounted  for  by  a  multiplicative  correction  factor 

Y  (zi )  =  1  -  exp  [(-jk1  V-q2_1  _  -±2  )  (h+  zY)  ]  (22) 

A  -  tt2 

i.  e.  , 

(£L)  +  E2(°'  (rL)  *  Y  (aj)  (r^  (23) 
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Fig.  2.  Ray  configuration  relevant  to  the  calculation  of  G2. 
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and 


Si  '  —i^  +  Sir  '  —i^  **  y  (zi>  Si  p 

.  (o) 


(24) 


where  G,  and  E'  are  given  by  Eqs.  (18,  19)  respectively.  Strictly  speak- 
ing,  neither  Eq.  (21)  nor  Eq.  (24)  are  correct  since  vertically  polarized 
field  constituents  do  not  experience  a  phase  inversion  upon  reflection.  These 
results,  however,  are  applicable  in  our  case  since  both  transmission  and 
reception  are  horizontally  polarized  and  depolarization  due  to  either  co¬ 
herent  or  incoherent  scattering  is  expected  to  be  small* 

In  the  analytical  sequence  which  follows,  ground  contributions  associated 
with  G00  (Figure  Id)  are  ignored  for  the  sake  of  convenience.  These  con- 
tributions  may  be  readily  accounted  for  by  a  procedure  identical  to  that  pre¬ 
sented  below.  However,  since  ground  effects  are  of  secondary  importance 
to  the  determination  of  the  spectrum,  as  is  demonstrated  below,  we  do  not 
return  to  investigate  this  contribution.  From  Eqs.  (18,  19,  20,  23  and  24) 
one  has 

[Si  <£’  £i>  +  •  S2  *  [^20)  {— t}  +  -20' 


2  z  y  (Zj)  y  (zT)  F(cp  ) 
8tt2  (T12-1)  pq4  ^  -r 


,  z  2 
exp  j  -jkj  [  |  P. |+  |P.-£TI+  f - V  T)2"1  (Z!+ZT)] 


T  1 
z,  +  z , 


-*2  lli-ITl+a  7==^  -  al-i  --Tl  } 


si : 


(25) 


where, 

R=  U 


Tl 


(1-llHl-ll)  n  r,  (£i  -  {li  -  ^  , 

—  J  *  U - : - T7 - ]•  £, 


(26) 


-1  -T 
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If  (a)  |  r  |  >>  |_r  ^  -  r_^  |  and  (b)  ignoring  vertical  (crossed  coupled)  components 
not  sensed  by  the  receiver  one  obtains  (see  Appendix  B  and  Fig.  3) 


•  2  .  .  -1  "  -T 


£  M  Sin  X  (cos  x  =  ‘ 


-1  "  -T  I 


(27) 


Furthermore,  in  the  transition  from  Eqs.  (18,  19)  to  Eq.  (25)  the  geometric- 
optical  path  (Lj)  has  been  replaced  by  the  range  (p  )  in  the  amplitude  factor 
(see  Fig.  3).  Similarly,  corresponding  to  the  second  term  in  the  integrand 
of  Eq.  (17), 

[Gj  (r,  rT)  +  G1R  (r,  r  T)  ]  •  CF,  (rT>  •  [E^o)  (r^  +  E. ^°]  ( )  ]  ~ 


A  2z  V  (Zj)  Y  (ZT)  F  (cp) 
E  8tt 2  cn2- 1)  pq4  lij  -iT 


(  z2 
exp  |  -jkj  [|£-£j  |+  |£-£tI  +  “-  -  (Z!  +  ZT)J 


Z1  +  ZT 

-Jk2  111  "  It  I  +  a  ,  '  a 

i  -  n 


111  It  I  } 


(28) 


A 


Under  conditions  (a)  and  (b)  above,  p  %  p  and  expressions  (25)  and  (28)  are 
identical.  The  substitution  of  Eqs.  (25,  28)  into  Eq.  (17)  together  with  the 
approximate  forms 


|p-PiI^£o-£o  '  £r  Lp-£tI  ~P0'£o*  ELt 


(29) 


results  in 


i  3  2  Z  Y  (zt)FM  /  z2  r 2 

^  — T  (l  1 1)  ~  -j%UQ  2  2  4  eXP  {-2Jki  P0-Jki  — +  jk!^  -1  ZT  + 

2tt  ( T|  -  1 )  p  o 


a 


v  1  -  T) 


-2 


+  Jk!  P_^  *  £T  \l  (1' 


t) 


(30) 


1 1 


18-4-  13471 


SCATTERING  POINT 


Fig.  3.  Top  view. 
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where 


3  2  If 

:(£_>  t)  =  J  d  r  £(£1(  t)  sin  \y  (z  )  - - ■  exp  <  jkj  •  £ 

V 


+  a 


N1 


.-2 


4  jkj  VT12-  1  zt  -  (jk2  +  a)  IjTj  -  £t|  | 

The  plane  wave  expansion 

-1  A  j 

exp  [-jk2  | _r  -  £T  |]  =  - j"  J  d 

2tt 


(31) 


Hi  -lT 


expf-jx.  Irj-^IJ  a 


2  T  2 

x  -  k2 


k2  =  k2  -  ja 


(32) 


will  be  found  convenient  in  the  forthcoming  analytical  sequence.  For  any 
plane  wave  constituent  (characterized  by  the  wave  number  ft)  the  angle  y  is 
constant  depending  exclusively  on  ft.  The  substitution  of  Eq.  (32)  into  (31) 
results  in 

3  exp  [jx  *  £  ] 

I  =  - —  J  d  x  - 2 - 7Tz -  sin^  X(K)  j  d  r ,  e  (r,  ,  t)  y  (z  ) 

2tt  x  -  k2 


exp  j  -j[  x  -  kj  -  kL  vV  -  1  zj  *  £j  ^  (33) 

Eq.  ( 30) ,  whenever  applicable,  constitutes  a  convenient  starting  point  for  a 
statistical  analysis  of  the  backs catter ed  wave.  Of  primary  interest  are  the 
correlation  properties  of  the  return  signal  (both  spatial  and  temporal),  the 
corresponding  spectra  and  the  mean  return  power.  The  next  section  is  dedi¬ 
cated  to  the  study  of  the  temporal  correlation  and  the  associated  spectrum. 

A  simple  and  readily  inter  pretable  result  is  derived. 

D.  THE  TEMPORAL  CORRELATION  COEFFICIENT:  C«(t) 

E 

The  field  temporal  correlation  coefficient  is  defined  by 


13 


I  —i  t(— ’  I  >  C£*T'  =  <§iT*-’t*  '  -IT  (£.’  t  "  T)>=  5i  t(— ,  t)  *  5lT  t_T^ 


-IT 


.  (o) 
-IT 


.(<>)*, 

-IT 


+  <  AE^  t(£,  t)  •  A  E1  T  (r_,  t-T  )  > 


(34) 


4az, 


Explicitly, 

,  >  lns*  ?  ?  2  7,4 \y (ztH 4 lF(cp>  lZ  — “T 

^(t)  •  <t-T)  =  «»0  HQ  |P|  2  2  2  5 - exp[j2kl£o.  VT  +  — -] 

n  ( T\  -1)  pQ  'l-T)"2 


(35) 


and 

-IT 


<AEirp(t)  •  AEf^  (t-T  )>  =  u,%„4  |p|2 


-IT 


o  o 


Z4  |y(zT)  |2F(cp)  \Z 

~  4  .2  ..2  8 

4tt  (T)  -1)  pQ 


<I(t)  I'  ( t- T ) > 


exP  [jk!  ‘  ZT  + 


2a 


n  -if2 


(36) 


where  the  relation 


£T(t)  -  £T(t  -  T)  =  v  t 


has  been  utilized.  Consequently,  one  has 

2  z4  |y(zT)|4|F(cp)|: 

<I§1TI  >  c,iTt  =  .‘ V  |e| 


o  o 


22  2  8 
n  (1)  -1)  Pq 


(37) 


1  + 


4  2 

(JO  |i 
o  o 


4tt2  |y  (ZT)  I2 


<I(t)  I  ( t  -  T )>  exp  [-jk^  ^  •  vt  -  2a 


A  -Tf  2 


]  exp  [j2k  p  •  vt  +  - ]  |  ( 38) 

l-°  -  / - —  J 

V  1  -T) 


with 

<I(t)  I*(t-T)>  = 


2  2 

1  Pi3  3  sin  X  (Hj)  sin  X  (k2) 


;d  1,1  Jd 


exp  jfKj  '  £p(t)  -  K2  •  £p(t-T)  ]  J  d3  r  J  d3  r2  C  (£  -  r_2  |  t)  y  (z^y  "(z^ 

V  v 

r-j -  a(z  +z  K 

exp  J-jKj-  +  jK2  -  £2  +  jki  ^  •  (r_i  -  r_2)  +  jk^ T)  -1  (Zj  -  z^)  j.  (39) 

V1  -  T)-  2 


14 


and 


Cg(£l  -  £2  I  T  )  =  <e  (£1 ,  t)  e  (£2,  t  -  T  )> 


r-  ,  ,  |  -in  •  (r,  -  r?)  d  x 
=  J  tfg  (H  I  7  )  e  J-  V-1 


(40) 


Eq.  (40)  is  based  on  the  assumption  of  statistical  homogeniety  and  stationar- 
ity  of  the  process  e  (r,  t).  The  following  change  of  variables  is  introduced 
for  convenience: 

A  1  +  1  A 

L  =  £i  -  Z2*  -  =  2  (£l  +  l2)  °r  -1  =  ^  -  2  - 

2 

^  1  .  ,  -  x  1  A 

—  ~  ~  —2‘  —  ~  1  (—1  +  —2}  or  —  1  “  —  “  ~2  — 

2 

£t  =  £T(t)  -£T(t  -  T )  =  V T  ,  £t  =  J  [£T(t)  +£T(t  -  T)  ] 


o  r 


V  T 


V  T 


£T(t  T)  £T "  2  -  £T(t)  I-T  +  2 


A  v  1  ,  „  +  1  A 

2  =  z  1  ”  2  2  ’  Z=  2  (Z1  +  Z2)  or  z1=  Z±  r  z 


(41) 


and  Eq.  (39)  becomes 

<I(t)  1  (t  -  t)>=  — L-  J  d  3  xJd^Jt 

4tt 


2  2 
sin  xT  sin  X2 

- 2 - - - 2 - “ — 

I H  +  i  H  |  -k2l  [  |h-±h|  -k-"2] 


exp 


A 

VT  4  j  K  ■ 

— 1 

d3R 

J 

V 

V 

A  A 

K  ■  r  -  K  « 

R  +  k1 

£0  • 

l 

,3  A  „  / A 
1  r  Ce(l 


1  A  * 

+  2  z)Z 

2  a  Z 

v  1  -  T)'2 


1 

2  Z) 


(42) 


here, 


1  A 

Xj  =  X  (H  +  2  *) 
2 


(43) 
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For  simplicity,  we  further  assume  that  y  does  not  vary  substantially  over  the 
correlation  length  (a),  i.  e. 


k,  a 


, : t~ 

Ti  - 


i  <<  i, 


a 


a  <<  1 


■Ti 


T| 


-2 


(44) 


so  that  the  approximate  relationship 

y  (Z  z)y  "  (Z  -  I  z)  «  \  y  (Z)  | 


(45) 


may  be  utilized  in  conjunction  with  (42).  To  the  same  order  we  may  also 

n — 

ignore  the  refractive  phase  term  k^  vT]  -  1  z  in  Eq.  (42).  It  is  noted,  how¬ 
ever,  that  this  approximation  is  unessential  and  may  be  relaxed  if  necessary 
With  (45)  in  mind  and  for  an  effective  scattering  volume  characterized  by  a 
linear  dimension  large  compared  with  the  correlation  length,  Eq.  (42)  yields 


<I(t)  I  (t-  T)  > 


(2rr y 
4rr4 


2 

sin  x 


.  2 


-  1  - 

—  ~Z  — 


-k22]L 


1  A  |  **2n 
k-2  k|  -k2  ] 


f4)  #e  (ki  Eo  +  k!  Jr}2  -  l  ~  1L  I  T  )  exp  [j  x  '  v  r  +  jv  1  £_  T  ]  (46) 

where , 

F(i)  =  J  d  3R  |  y  (Z)  |  2e  ^  (47) 

yi-r)-4 

and  is  the  spatial  spectral  density  defined  in  Eq.  (40).  Eq.  (46)  may  be 
cast  into  the  form 

<  I  (t)  I*  (t- t)  >  =  ~  exp  [  j(k^  +  k^  Jr}2- 1  z^)  •  v  r  ]  J  d 2  v  F(v)M(£_  |r) 

exp  [  j  C_.  r_  T  ]  (48) 
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where 


M(v_  |T  )  =  d  v 


2  2 

0£  (h  T  )  s in  X1  sin  x  2  exP  t -j  H_  •  VT  ] 

1  .  "|2  7  2-,  r  I,  1  „  i2 


C  lk2P0  "  *  +  J  1  [  lk2£-o  ~k  —  ^"k2  (49> 


Xt  X(k2^-i±^l) 
Z 


and  a  new  variable 

*  kj  e^,  +  kj  .4? 


(50) 


-  1  z  -  H  k 
— 1 o  — 


2^>  -  SL 


(51) 


v4^i"  £ 


has  been  replaced 


has  been  introduced.  The  parameter  + 

2  2  21/2 

by  its  proper  magnitude  (k^  =  kj  +  k^  ( 1  -7]  )  ]  )  and  approximate  direction. 

For  sufficiently  large  volumes  (this  time  compared  to  the  wavelength) 
and  for  slight  to  moderate  losses  F(h)  is  sharply  peaked  at  and  Eq.  (48) 

reduces  to 

■a  <3  •  a 

<r  I(t)  I*  (t  -T  )  >  as  ~  exp  [jk^E-o  *  v_  t  1  C  J  cl  h  F(h  )  e  1  —  —  T]M(0|t)  (52) 

where  the  bracketed  term  is  readily  evaluated  as 


r.3 


d  v  F ( v)  e 


J  v  .  L 


T 


u  *  f  d3  R  I  Y  (Z)  |2  exp  ]6  (R  -  r_  )  = 

V  yi-Ti  T 


2  az 


T 


( 2 tt  )  |v  (z  )  I  exp  [— - T 

yi-T) 


(53) 


The  substitution  of  Eq.  (53)  into  (52)  and  subsequently  into  Eq.(38)#  yields 
1  +  4uu  4p  2M(0  | t  ) 

Ce(T|“  ,  a  °«  n,  eXP  C2jk'  '  -T]  <54) 

1  +  4a;  |j,  M(0t=0) 

o  o 
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with 


M(0  |t)  =  j  d3n 


3  0£(h  |t)  sin  x  (k2^0Z^)  exp  [ -jn_  •  v  t] 

"  |k2  ^-o  "  —  I  ”  k2  J  "  lk2  ^-o  ”  —  I  ~k2  - 


(55) 


v_  is  presumed  horizontal  as  has  been  previously  stated.  It  may  be  observed 

that  since  0^(y  ,7)  is  symmetrical  in  r  it  follows  that  M(0  |  7)  =  M(0  |-f)  and 

the  anticipated  constraint  C^(t)  -  C*  (-7)  is  satisfied.  It  is  furthermore 

E  E 

2 

observed  that  upon  the  cancelation  of  the  factor  y  ( z  ^ )  in  Eq.  (53)  with 
that  appearing  in  Eq.  (38),  ground  effects  totally  disappear  to  the  cited  order 
of  approximation.  The  lack  of  spectral  sensitivity  to  the  presence  of  the 
ground  -  vegetation  interface  is  traceable  to  the  fact  that  the  targets  motion  is 
horizontal  (parallel  to  the  interface).  Deviation  from  this  configuration  may 
result  in  a  substantial  increase  of  sensitivity  to  ground  effects. 

We  now  attempt  the  evaluation  of  the  integral  in  (55)  subject  to  the 
simplifying  assumptions: 

(a)  v  =  v  (motion  toward  or  away  from  the  receiver). 

(b)  0  (h  |t)  =  0  (k  |t),  h  =  |  h_  1  (statistical  isotropy). 

The  integration  process  is  carried  in  the  spherical  coordinate  system  de¬ 
picted  in  Fig. 4.  One  has, 


k 


2^0  -5- 


4-  h 


2k^  h  |i  ,  g 


cos  0 


or 


'k2  ^  "  * 


=  (H  -  Hj)  (K  -  H2) 


(56) 


(57) 
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4 


Fig.  4.  Spherical  coordinate  system  in  wavenumber  space. 
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whe  re 


[  1  +  J\  -  — 


H1  -  k2*J'  tl+  a/1  -  2  ] 

2  k2  p 


(58) 


so  that  (55)  becomes 


1  “  ?  0  (h  |t)  sin  x(M- »9)  exP  [  -jn  v  (i  t  ] 

M(0  |t)  =  dcp  j  dp,  J  dHH  - * - * — 

0  -1  0  (k  -  Hj)  (h  -  H^)  (h  -Hj  )  (k  -  h2  ) 


(59) 


The  integrand  in  (59)  has  the  property  f(n ,  p )  =  f(-H,  -p  )  from  which  it  follows 
that 

-  1  OO  1  —00 

j  dp  J  dn  f  (h  ,  p  )  =  d  p  dn  f  (H  ,  p  ) 

1  0  6  00 

The  limits  of  integration  in  Eq.  (59)  may  be  changed  accordingly  and  the 
H  -integration  is  readily  performed  by  a  residue  calculation.  The  correspond¬ 
ing  deformation  of  the  integration  path  is  depicted  in  Fig.  5  and  leads  to  the 
result 


2tt  1^ 

M(0  |t  )  =  2nj  J  ^  |  dp  sin  \(m<»cP) 
0  0 


#  2  „  -,v 

H1  ^H1  lT  )  exP  C  -JM-  t  VpT  ] 

(Hj  -h  x  )  (h  l  -  h  2 )  (K*  -*2*) 


k2  ^e^Z  lT  ^  exp  ^  "J  h2  V(j,t  ^ 

ft 

(H 2 -H  1  >  (H 2  ~ H  1  )  (H2  ”H2  ^ 


(60) 
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18-  4-13473 


Fig.  5.  The  k -plane  description. 
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if  t  >  0  and  to 

2tt  1 

M(0  |t)  =  2rr  j  dcp  |  dp,  sin  x(M->9) 

0  0 


2 

0£(h1  |t  )  exp  [  -  j  H  j  v(i  t  ] 
<H  1  ~H2>  (H  1  ~H  1  *  ^H'l  "H2  ^ 


■><•2  *  | 

H2  «Pe(H2  ^  GXP  " j H  2  VM,T  ' 

-  —  5f  £ 

<w2  “Hl)  (h 2  “H2)  (h 2  "H1  } 


(61) 


if  T  <  0. 

The  substitution  of  the  explicit  forms  of  h  ,  given  in  Eq.  (50)  leads  to  rather 

i  ,  L* 

cumbersome  expressions.  Although  the  analysis  can  be  pursued  further  with 
these  expressions  we  concentrate  on  the  limit  of  small  but  nonvanishing  loss 
(per  wavelength)  in  which  the  following  simplified  forms 


► 


H1  ~  j  f  ’  H2  =  2k2  f-  (62) 

are  obtained  retaining  first  order  terms  in  —  .  Obviously  the  expansion  of 
^  is  near  (j,  0  where  the  exact  forms  should  be  used.  How¬ 

ever,  this  corresponds  to  a  small  fraction  of  the  range  0  <  g  ^  1  which, 
in  turn,  is  associated  with  a  small  contribution  to  the  ensuing  integral.  The 
contributions  associated  with  the  residues  at  y  *  and  y^  (first  terms  in 
Eqs.  (60)  and  (6l)  respectively)  are  proportional  to  a  and  may  be  neglected. 
The  only  cp  -  dependent  factor  in  Eqs .  (60,6l)  is  (see  Fig.  4). 

sin^X  (m*  >9)  =  [  1  -  sin2  cp  sin2  6]  2  =  [ 1 -( 1  -  g  2)  sin2cp]  2  (63) 

which  is  readily  integrated  over  cp ,  yielding 


2tt  4 

£(M- )  =  f  sin  X(M>  »<P)  dcP 
0 


J  [1-(1-^2)  sin2  cp  ]  2  dcp  =  j-  [  3  +  2g2  +3p4]  (64) 
0 
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Retaining  loss  effects  to  lowest  order  and  ignoring  regions  near  g  ^  0  so  that 


Z  - 1 

a  «  2k0g  and  a  «  (vt  ) 

Z  c 

one  obtains  ,  ~ 

1  -j2k  [i  v  t 

M(0  |t)  —  dp,  p  /3(n)  0£  (2k  h  |t  )  e  (65) 

0 

We  observe  that  Eq.  (65)  is  not  explicitly  dependent  on  the  scattering  volume 
(V)  as  would  be  the  case  in  a  lossless  case.  This  is  indicative  of  the  fact 
that  the  effective  scattering  volume  is  basically  determined  by  the  losses 
(Appendix  A).  The  last  observation  does  not  contradict  the  low  loss  assump¬ 
tion  introduced  previously,  a  may  indeed  be  small  per  wavelength  or  com¬ 
pared  to  (v  t  )  ^  and  at  the  same  time  be  sufficiently  large  to  effectively 

restrict  the  scattering  volume  (V  in  Fig.  A.  1  is  typically  very  large). 

8 

Let  us,  for  the  time  being,  ignore  the  motion  of  the  scattering  centers 

represented  by  the  t  dependence  of  0^.  We  return  to  study  this  effect  in 

Section  E,  below. 

2 

Let  u  g  denote  a  new  integration  variable.  Eq.  (65)  may  be  cast  into 
the  form 

1  _ 

M(0  |t)  =  ~  J  d  u  j8  (  J u  )  0£(2k2  Ju  )  exp  [  -2  jk2  uv  t  ]  (66) 

and  the  Fourier  inversion  of  Eq.  (54)  results  in 
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(67) 


or 


(68) 


for  0  <  a)  <  2k^v  and  equals  to  zero  in  the  complementary  part  of  the  spec¬ 


trum.  The  superscript  (o),  above  indicates  the  fact  that  the  motion  of  the 
scatters  has  been  omitted.  We  terminate  this  section  on  the  following  notes: 


(a)  The  singular  term  5(cu-2k^  v),  representing  the  direct  target 


return,  will  be  generally  diffused  owing  to  such  factors  as  ran¬ 
dom  or  systematic  variations  of  the  target  velocity,  finite  time 
broadening  etc.  It  basically  represents  the  target  spectrum  in 
the  clear  (i.  e.  ,  in  the  total  absence  of  multipath  effects). 

(b)  The  second  term  represents  single  scatter,  multipath  effects. 

It  increases  with  decreasing  a  and  disappears  as  a  increases. 
It  could  be  anticipated  that  the  spectral  return  from  dense, 
highly  lossy  vegetation  slab  will  resemble  that  obtained  in  the 
clea  r . 
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(c)  As  a  decreases,  the  relative  significance  of  multipath  effects 
increases  to  the  point  in  which  the  effective  scattering  volume 
is  no  longer  loss  limited.  This  limit  is  irrelevant  in  the  pres¬ 


ently  discussed  physical  context  but  is  of  great  interest  in 
most  other  situations. 

(d)  Eq.  (68)  constitutes  a  simple  and  readily  interpretable  result. 
The  interpretation  is  presented  in  Section  F.  It  facilitates  the 
prediction  of  0  ^°'  (uu )  provided  that  the  spatial  spectral  density 

iij 

function  (0^)  is  known.  Conversely,  Eq.  (68)  may,  under  favor¬ 
able  circumstances,  form  a  basis  for  the  predictions  of  0^ 

from  the  measurement  of  the  target  spectrum  (0  ^°^)  under 

E 

no  wind  conditions. 


E.  EFFECTS  DUE  TO  MOTION  OF  THE  SCATTERING  CENTERS. 

The  question  of  spectral  perturbation  caused  by  the  random  motion  of 
the  scatterers,  characterized  by  a  velocity  field  U_  (t),  is  discussed  in 
Ref.  1  where  the  properties  of  the  clutter  return  are  investigated.  The  anal¬ 
ysis  is  based  on  two  major  assumptions:  (a)  the  random  processes  £(r_) 
(describing  the  spatial  distribution  of  the  scatterers)  and  U  (t)  are  statisti¬ 
cally  independent  and  (b)  e(r^,  t)  is  a  conservative  process,  i.  e. 

t  +  T 

e(r_,  t  +  t)  .  e(r_-  J  U  (s)  ds,  t)  (69) 

t 
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Consequently,  the  space-time  correlation  function  defined  in  Eq.  (40)  is 
describable  by 


C 


E(-l 


(70) 


where  C  is  the  spatial  correlation  and  <  denotes  averaging  over  the 
velocity  field.  The  spatial  Fourier  transformation  of  Eq.  (70)  results  in, 


t  +  t 


0e(K_lT)  =  - }~  J  d3  r  e'J~  ’  -  <C£(£  +  j  U(s)ds)>u 


(2tt)' 

1 

(  2tt  )  ^  £  t 


.  I  ~ 


t  +T 


d  r'  e^—  —  C  (r_! )  <exp  [ -j  h_  •  j  U(s)ds]> 


=  (H_)  W  (h_|t) 


(71) 


where 


t  +T 


w(hJt)  =  <exp  [  -  j  h_  *  U(s)ds]>u 


(72) 


t  +T 


,  r_  =  r^  -  _r  r_l  =  r_  +  J  U(s)  ds  and  0  is  the  Fourier  transform  of  C 
1  1  t  ~  6  6 

i.  e.  the  spatial  spectrum  in  the  total  absence  of  motion.  The  substitution  of 

Eq.  (71)  into  Eq.  (55)  leads  to 


4 

0  (h)  sin  Xw(h_1t)  exp  [  -j_H  •  v  t] 


M(0  |t  )  =  J  d  h 


3  ^0 


-  k/]  [  |k2£^-H.j2-fc^2] 


(73) 
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If  the  random  velocity  field  U  (t)  is  isotropic  it  follows  that  w  is  exclusively 
a  function  of  h  =  |  H_]  .  It  follows  from  Eq.  (73),  via  a  procedure  analogous 
to  that  leading  to  Eq.  (65),  that 

1  ^ 

M(0|t)  =  —  J  d)i  )i  i3()i  )  0e(2k2M, )  w(Zk2i^  I  t  )  exp  [ -jZk2M' "  v  t  ]  (74) 

or  analogously  to  Eq.  (66) 

1  _  _  _ 

M(0|t)  2^-  J  du/3(y'u)  0£  (Zk2  Vu  I  t  )  w(2k2Vu  |  t  )  exp  [ -jZk2uvr]  (75) 

The  Fourier  transformation  of  Eq.  (58)  becomes, 

0E(w)  [  1  +4  wo4po2  M(0|t=  0)]'1  {6(cu-2k1v)  +  4u)o4goZ  £  J  du  0(*/u  ) 

o 

0^(Zk2  vu  )  W(Zk2  */u  |cu-Zk^  v  +  Zk2  v  u)  j-  (76) 

where 

—  i  +  00  _ 

W(  Zk2  J\\  |  cu )  =  2^—  >  dT  w(Zko  v  u  |  t  )  exp  [  -jcuT  ]  (77) 

_  00 

Eq.  (77)  reduces  to  (67)  as  W(2k7  a;u  |<ju)  approaches  6(uu). 

The  nature  of  W(  2k  ^  Vu  |  )  has  been  investigated  in  Ref .  (  1 )  whe  re  two 

alternative  models  have  been  suggested.  The  first  of  which  involves  the  as- 
t 

sumption  that  U  (s)  ds  constitutes  a  multivariate  Gaussian  process  pos- 
o 

sessing  the  features: 

(a)  <  U  >  =  0. 

(b)  Orthogonal  velocity  components  are  uncorrelated  i.  e. 

<  U.  U.  >  =  0,  i  4  j. 
i  J 

2  2 

(c)  The  velocity  field  is  isotropic.  Specifically ,  <U  >-<l > 

1  2 


2 


One  obtains,  via  Eq.(72), 


2  2  7 

w(H  |t  )  =  exp  [  -  ct  H  f  (T-S^CU  (s)  ds  - 
u  0 

where 

C  (s)  =  3  a  "2  <  U.  (t)  U.  (t  +  s )> 

U  U  11 


(78) 


(79) 


is  the  correlation  coefficient  which  is  independent  of  i  owing  to  the  presumed 
isotropy  of  U  (t).  Qualitative  and  quantitative  details  concerning  the  form  (78) 
of  w(h  |t  )  are  presented  in  Ref.  (1),  so  that  further  elaboration  seems  unnec- 
sary.  An  alternative  model  (Ref.  1)  views  the  velocity  field  as  quasiharmonic, 

i.  e. 

U  ( r_,  t)  =  £?  (rj  (r)  cos  [f?(r_)  t  +  cp(r)  ]  (80) 

where  and  cp  are,  assumedly,  mutually  independent  random  vari¬ 

ables.  The  substitution  of  (80)  in  Eq.  (72)  and  upon  Fourier  transforming 
the  result  one  obtains 


W (h_  |uj )  =  S  <J  (h_  •  0  )  6  (U>  -  nfi)  exp  [ -jncp  +  jn_  •  psincp]>  ~ 

n=  -  30  r  >  > 


(81) 


where  the  expansion 


exp  [  -j  H_  ’  j>  sin  (fiT  +cp  )  ]  =  E  Jn(n_  • 

n=  -» 


-jn  (Qt  +  cp) 


(82) 


has  been  utilized  and  <  >  n  denotes  sequential  averaging  over  the 

random  variables  cp  ,  Cl  and  respectively.  Assume  the  cp  is  uniformly 

distributed  in  the  interval  ( -  TT  ,  tt  )  and  let  P(fJ)  denote  the  density  function  cor¬ 
responding  to  Q .  Furthermore,  assume  |  (the  maximum  value  of  the 
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projection  of  £_  along  h)  to  be  Rayleigh  distributed  with  the  locations  of 

its  peak.  One  obtains  (see  Ref.  1) 

00 

exp[HZpm2]W(H  |u>)  =  ■  Iq(h2P  f)6(«u)  +  Z  In(H2pm2)  [P(^)+P(-|)]  (83) 

n  -  1 

where  I  (x)  are  the  modified  Bessel  functions  of  the  first  kind.  The  follow- 
n 

ing  functional  dependence  for  P  (fi)  has  been  found  suitable, 

(Q/Q  )p 

P(0)a  - 2 - -  (84) 

i  +b(o/n  )q 

o 

where  Q  ,  b,  p  and  q  may  be  determined  experimentally  from  the  spectrum 
of  clutter  return  at  low  winds. 

F.  THE  TARGET  SPECTRUM:  HEURISTIC  CONSTRUCTION. 

The  Doppler  shift  associated  with  the  scattering  events  depicted  in 
Fig.  6  are,  neglecting  relativistic  effects,  given  by 

uu  =  k  v  (cos  j9  -  cos  a  )  (85) 

where  uu ,  k,  and  v  denote,  respectively,  the  Doppler  shift,  the  wavenumber 
of  the  incident  radiation,  the  target  velocity  (presumed  constant),  a  and  /3 
are  the  angles  between  v_  and  the  wavevectors  characterizing  the  waves  in¬ 
cident  upon  and  scattered  by  the  target,  respectively.  Fig.  6a  depicts  a 
configuration  in  which  the  wave  undergoes  a  random  scattering  event  having 
been  scattered  by  the  target,  while  Fig.  6b  describes  the  converse  situation 
in  which  the  random  event  occurs  prior  to  the  target  scattering.  The  two 
cases  result  in  identical  shifts  given  by 
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(a)  a  *  it  ,  (3  -  -9t 


(b)  a  »  w  +  0,  ,  /3=o 


Fig.  6.  Doppler  shift  by  scattering. 
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uu 


(86) 


=  k  v [  l  +  cos  9s  ] 

where  0g  ,  the  s cattering  angle ,  equals  -0  in  Fig.  6a  and  to  a-n  in  Fig.  6b. 

It  is  noted  that  u)  is  never  negative.  Reversal  of  the  shift  direction  is  ob¬ 
tained  only  by  the  reversal  of  v.  This  is  a  direct  consequence  of  the  fact  that 
v  is  presumedly  directed  toward  the  transmitter  (receiver).  The  spectrum 
of  the  target  return  is  determined  by  the  mean  power  return  from  the  entire 
range  of  the  scattering  angle  0  .  Consider  a  scattering  volume  on  which  a 

plane  is  incident.  The  effective  scattering  cross  section  of  the  volume  is 

(2) 

proportional  to  v 

?  ~  0 

sin  X  0£(2k  sin  ~Y~  )  (87) 

subject  to  the  following  constraints;  (a)  single  scatter;  (b)  linear  dimensions 

of  V  large  compared  to  the  correlation  distance  characterizing  G(r_)  in  V; 

(c)  e(r_)  is  statistically  homogeneous;  (d)  e(r)  is  statistically  isotropic; 

.  ||  -  .  Z 

(e)  "fa  r  field"  interaction.  An  additional  "pola r ization  factor  sin  \  appears 
due  to  the  target  scattering  itself  and  (87)  becomes 

4  .  & 

sin  x  0£  (2k  sin  ~r  )  (88) 

The  angle  X  is  n°t  exclusive  function  of  0  .  For  a  fixed  value  of  0  ,  the 

s  s 

4 

factor  sin  x  may  be  "averaged"  over  the  angle  cp  depicted  in  Fig.  4,  since 
variations  in  cp  (with  0^  constant)  do  not  correspond  to  spectral  variations. 
The  result  of  this  averaging  is  the  function  /3/2tt  defined  in  Eq.  (64).  The 
substitution  of  (from  (86)) 

sin  "T"  =  '  Tfe  (89) 
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leads  to  the  expression 


0(i- 


uu 

2kv 


)  0e(2k  jr~- 


SL.  ) 

2kv 


(90) 


which  is  identical  (disregarding  proportionality  constants)  to  the  "multipath* 1 

term  appearing  in  Eq.  (68). 

G.  THE  MEAN  POWER  RETURN 

While  the  main  thrust  has  been  directed  toward  a  gain  of  insight  into 

the  mechanisms  which  control  the  spectral  characteristics  of  the  target 

return  an  estimate  of  the  mean  power  return  is  of  great  interest  as  well. 

The  desired  results  stem  directly  from  the  non -normalized  forms  (35-36) 

evaluated  at  t=  0.  For  a  matched  antenna, 

2  2 

<pt>  =  bf g(:p)  <s>  "  br  g(cp)  £ 


g2(y ) 

p  ,-> 

o  (2n ) 


list 


(y(zT)  I 


(ti2-D2po8 


■exp  [ 


4azr 


A  -if 


—  ]  C  1  +  4ou  4p  2  M(0  |0)  ] 
2  °  o 


(91) 


where  P  ,  g  (cp )  and  <S>  denote  the  transmitted  power,  the  antenna  gain  at 
0  0  and  the  Poynting  vector  at  the  receiver,  respectively.  M(0  |0)  and 

y(z  )  are  given  in  Eqs.  (66)  and  (22),  respectively,  and  jS(y^u  )  in  (66)  should 
not  be  confused  with  the  constant  polarizability  factor  introduced  in  Eq.  (1) 
and  appearing  in  Eq.  (91),  above.  The  first  and  second  terms  in  (91)  are 
recognized,  respectively,  as  the  contributions  associated  with  the  direct 
and  multipath  returns.  The  sensitivity  of  the  mean  return  to  medium 
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(scattering  and  absorption)  losses  (represented  by  the  factor  exp  [ 


4a  z , 


2  y  i  -Ti 

to  ground  proximity  (represented  by  |  y  (z^  j  ),  to  range  (represented  by 
-  8 

p  )  and  to  antenna  height  above  the  fluctuating  slab  (represented  by  the 
4 

factor  z  )  are  self-evident. 

The  signal  to  clutter  power  ratio  is  determined  via  a  comparison  of 
Eq .(91)  with  1(127)*  resulting  in 


]). 


P  >/<"P  > 

T  clutter 


10  |V(<?)  |y(»T)  |4exp[2a(2zT+h)]  [  1  +4wo\2M|0|  0)]  (g2) 


v0dh[f^]d?]  [sin  ‘0e(2kr2klyr 


i) 


In  (92)  h  and  d  denote  the  effective  forest-slab  height  and  the  depth  of  the 
range  cell,  respectively.  The  first  and  second  arguments  in  0£,  above, 
refer  to  the  respective  horizontal  (  |jjt^  |  =  |  H  -  z^  h  |  =  2k  ^ )  and  vertical 


(h  2k 
z 


./ 


-1  )  components  of  h  . 


- - -  2 

*Eq  .1  (127)  c ontain  s  two  typographical  errors:  (a)  the  factor  (T]  -  1 )  in  the 

denominator  should  be  squared  and  (b)  the  sin  h  argument  should  read 

2ah/  /l-Tf2  . 

v 
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APPENDIX  A: 


THE  EFFECTIVE  SCATTERING  VOLUME 


The  volume  which  actively  participates  in  the  scattering  process 
described  in  the  text  is  defined  by  the  intersection  of  three  factors: 

(a)  the  finite  extent  of  the  slab,  (b)  the  radar  gating  and  (c)  the  background 
(effective)  loss,  (a)  above  requires  no  further  comment.  The  effect  of  the  timing 
constraints  (b)  is  depicted  in  Fig.  A.  1  in  which  ground  effects  are  omitted 
(but  could  be  readily  accounted  for  by  appropriate  imaging).  The  scattering 
volume  (V  ),  defined  consistently  with  the  radar  gating,  constitutes  a  para- 
boloid  of  revolution  (with  zn  the  axis  of  symmetry)  determined  below.  Con¬ 
sider  the  cross  section  of  the  scattering  volume  with  the  plane  yn  0 
(Fig.  A.  1).  The  phase  accumulated  by  the  wave  incident  on  the  target  having 
undergone  a  scattering  event  at  the  point  (xn,  z")  is  a  function  of  (xn,  zM) 
given  by 


k2[  Jx 


-  z' 


]  =  constant 


(A.  i) 


stemming  from  the  fact  that  the  incident  phase  front  (prior  to  scattering)  or 
the  phase  front  of  the  backs cattered  signal  (having  undergone  a  scattering 
event)  are  approximately  planar  (i.  e.  G^  is  a  constant  angle,  insensitive  to 
variations  in  0^).  Eq.  (A.  1)  yields  the  family  of  parabolas  given  by 


x"2  -  2cz"  =  C2,  c  >  0  (A. 2) 

or  in  the  x1  -  z’  coordinate  system: 

x"  ~  x’  cos  0  ^  +  z1  sin  0^,  z"  z*  cos  0^  -  x1  sin0^  (A.  3) 
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PARALLEL  TO  Z “  AXIS 


Fig.  A.  1.  Gating  limited  scattering  volume. 
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one  has 


( «/l  -T|  * 1  +  ^  ) 2  -  2  c  ( /l-Tl  2  z 1  -  =  c 2  ( A .  4 ) 

in  which  the  approximate  form  sin  0  0  7]  ^  has  been  utilized.  The  maxi¬ 

mum  value  that  c  may  get  is  determined  by  the  size  of  the  range  cell.  It  is 
not  hard  to  show  that  all  scattering  points  external  to  the  paraboloid  char¬ 
acterized  by  c  >  c  -  y  (d  -  the  range  cell  depth)  will  not  reach  the 

max  c, 

receiver  on  time. 

The  third  constraint  on  the  scattering  volume,  the  background  loss,  is 
probably  the  most  important  under  the  prevailing  circumstances.  The  loss 
limited  volume  may  be  defined  as  the  locus  of  all  scattering  points  (x1,  y1,  z1) 
obeying  the  inequality 

a  [  -L  +  L'  +  L"  ]  £  N  (A.  5) 


where  N  is  a  constant  representing  the  maximum  loss  difference  between 
the  paths  Ln  -f  L1  and  L,  (Fig.  A.  2)  beyond  which  we  may  totally  neglect  the 
scattered  wave  contribution,  a  is  the  loss  rate  defined  in  the  text  and  the 
geometrical  length  segments  are  shown  in  Fig.  (A.  2).  The  volume  bounds 
are  given  by 


ft  - 


cos  0 . 


+  [p'  + 


1 2  ,  i  2-, ^ 


:os  0 


-  =  N 
2  a 


(A.  6) 


or 


>,2  =  <  —  ♦ 


a  cos  6. 


i2  -  *.2  =  ,N 


■  wN  x  z* 


z'  =  (-  +  — -  z')(-  + — £-5-  +z')  (A. 7) 
va  cos  0 'a  cos  0O  '  '  ' 
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Fig.  A.  2.  Loss  constraints. 
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GROUND 


Fig.  A.  3.  The  effective  scattering  volume. 
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Approximately,  presuming  cos  0^ 


[l-Tf2]^  <  < 


one  has 


P 


cos  0  . 


(A.  8) 


The  effective  scattering  volume  is  depicted  in  Fig.  (A.  3)  in  which  the  angle 
3  is  given  by 

tan  6  =  757^7  «  [i-tT2]'4  (a.9) 

and  from  j  T)  -  1  |  <  <  1  it  follows  that  3  0  In  Fig.  (A.  3),  we  are  reminded 

again,  that  ground  effects  have  been  ignored  but  may  be  readily  accounted 
for  by  suitable  imaging  and  an  analogous  construction  sequence.  This  exten¬ 
sion  is  generally  important  for  ground  targets  but  is  insignificant  for 
spectral  considerations  as  long  as  the  target  motion  is  horizontal. 
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APPENDIX  B: 


THE  DETERMINATION  OF  EQ.  (27) 


(r  -  r  )  (r  -  r  ) 


lE  “Li 


<Ej  -Et)  (Et  -EtJ 

)  •  (1 - L - i - 75 - L  '  *o  =  (B-  1} 


1  "  — T 


cp  - 

■‘-O 


(E"E1)(E"E1)  •  E0  (E !  "Et) 


IE  “  E  j 


IE  i  "  Et 


cos  X  ^ 


~ — 1>^— ~ — i  >  *  (— i  -v 


■cos  x 


I—  “  —  1  *  ■  —  l  "  — T 


where  cos  X  ~  <£_o  *  (r^  ^  T  ^  I  L  i  ~  I* 

Under  the  assumption  _r  >>|r^-r^|we  have  (r_  -  r^  )  *  (£^%0and 
hence, 


^T7 


r  -  r 
-I  -T 


(r-r  )(r  -  r  )  r  -  r 

P“*-o  -  Cl 7-,r1'  ]  '  Ir-Tjl  "fr  ~r  |  cos* 

|_r-r,  |  —  — T  1  '—1  — T 


-1   -T 


cp  -  z  z  +cpcp  1  •  -| - r  cos  X 

-0-0  *o*-oJ  |r>1-ETl 


(B.  2) 


Ignoring  the  vertical  component  (not  sensed  by  the  receiving  antenna)  results 
in 

2 

R  es  c£_o  sin  X  (B.  3) 

Similarly,  under  the  assumptions  above, 


(E-ET)<E-Er)  n 
E"  Cl - -  T  ~  >  •  (~ 


.  (E,  -  ET)  ( LX  -ET)  rp 

- i - 2 - 1  9  1  -  )  •  ^o  (B.4) 


iE  -  Et 


lEi  Et 
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since 


(r  -  £_t )  (r  -  £t) 


^-o  2-o 


( B.  5 ) 


we  have, 


£  *  2-0  "  {l  "  £-o^o)  ‘ 


cos  X  =»  E  (B-  6) 
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